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Time-dependent  queues"  h 


ioseph  5.  Keller,  HAM  Rev., 


submi tted. 

A  single  server  queue  is  considered  having  exponentially  distributed 
inter-arrival  and  service  times  with  slowly  changing  *ime-dependent  rates 
A(ct)  and  u(ct).  The  parameter  is  the  ratio  of  an  average  inter¬ 
arrival  time  to  the  time  over  which  the  rates  change  appreciably,  so  it 
is  small.  Therefore  an  asymptotic  solution,  valid  ro r  small,  is  con¬ 

structed  for  the  time-dependent  queue  length  probability  distribution. 

It  consists  of  five  typical  parts  corresponding  to  five  typical  time 
periods.  They  are  the  initial  period,  the  period  of  light  traffic  when 
A/y  <  1,  the  saturation  transition  period  when  A/u  Increases  through 
unity,  the  oversaturation  period  when  A/y  starts  out  greater  than  unity 
and  then  decreases  below  unity,  and  the  transition  period  at  the  end  of 
oversaturation,  when  the  queue  returns  to  the  light  traffic  condition. 

By  combining  the  solutions  for  these  five  intervals,  the  solution  for  any 
queue  with  slowly  varying  rates  can  be  obtained,  borne  of  these  parts 
were  found  previously  by  G.  F.  Newell,  and  some  of  the  formal  results 
have  been  shown  to  be  asymptotic  by  W.  A.  Massey. 

63.  "Shock  profile  solutions  of  the  Boltzmann  equation"  by 
R.  E.  Caflisch  and  B.  Nicolaenko,  Comm.  Math.  Phys.,  submitted. 

Shock  waves  in  gas  dynamics  can  be  described  by  the  Euler  Navier- 
Stokes,  or  8oltzmann  equations.  We  prove  the  existence  of  shock  profile 
solutions  of  the  Boltzmann  equation  for  shocks  which  are  weak.  The  shock 


is  written  as  a  truncated  expansion  in  powers  o:  toe  shock  strcrutn,  toe 
first  two  terms  of  wm'cn  come  exactly  from  the  Taylor  tanh(x)  c'ofi le 
for  the  Navier-Stokes  solution.  The  full  solution  is  found  by  a  projec¬ 
tion  method  like  the  Lyapunov-Schmidt  method  as  a  bifurcation  from  the 
constant  state  in  which  the  bifurcation  parameter  is  the  difference 
between  the  speed  of  sound  Cg  and  the  shock  speed  s. 


64.  "Arrow’s 


-  >  ■> 


optimal  pricing,  use  and  exploration  of  un¬ 


certain  natural  resources"  by  P.  S.  Hagan,  R.  E.  Caflisch,  and 
J.  B.  Keller,  Econometrics,  submitted. 


Arrow's  model  for  finding  the  optimal  rates  of  exploration,  consump¬ 
tion  and  pricing  of  a  randomly  distributed  natural  resource  is  analyzed. 

At  first  the  model  is  modified  so  that  each  discovery  reveals  an 
arbitrary  amount  of  resource.  Then  it  is  analyzed  asymptotically  when 
this  amount  is  small,  and  analyzed  approximately  when  the  amount  is 
medium  or  large.  In  both  cases  explicit  results  for  the  optimal  explora¬ 
tion,  consumption,  and  pricing  policies  are  obtained.  Finally  a  numerical 
method  is  devised  which  is  applicable  whatever  the  size  of  the  dis¬ 
coveries.  The  results  using  it  are  found  to  agree  with  the  analytical 
results  in  both  cases. 
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65.  "Radiation  transport  in  a  hot  plasma"  by  Russel 


‘  :'l  i  sch  , 
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SIAM  J.  Appl.  Math.,  submitted. 

In  various  astronomical  events,  Compton  scattering  is  the  dominant 
interaction  mechanism  between  the  radiation  and  electron  fields  in  a 
plasma.  As  a  result  the  radiation  distribution  rapidly  becomes  a  Bose- 
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Einstein  given  by  A  (cie  -  1)  with  A  the  wavelength  and  a  >_  1,  a 
time  dependent  constant.  However  due  to  the  divergence  of  the 
Bremstrahlun.g  emission  spectrum  at  long  wavelengths,  the  distribution 
must  be  Planckian  there.  In  between  is  a  transition  region.  Based  on 
these  principles  an  expansion  is  found  which  is  an  asymptotically  exact 
solution  of  the  Fokker-Planck  equation.  Then  a  simple  equation  for  the 
time  evolution  of  the  radiation  temperature  is  obtained.  It  is  solved 
numerically  and  compares  well  with  previous  methods. 


66.  "Jets  rising  and  falling  under  gravity"  by  J.  B.  Keller  and 
J.-M.  Vanden-Broeck,  J.  Fluid  Mech. ,  submitted. 
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Steady  two  dimensional  jets  of  inviscid  incompressible  fluid,  rising 
and  falling  under  gravity,  are  calculated  numerically.  The  shape  of  each 
jet  depends  upon  a  single  parameter,  the  Froude  number  A  =  qc(Qg)  ^ , 
which  ranges  from  zero  to  infinity.  Here  qc  is  the  velocity  at  the  crest 
of  the  jet,  i.e.  the  highest  point  of  the  upper  surface,  Q  is  the  flux  in 
the  jet,  and  g  is  the  acceleration  of  gravity.  For  A  =  °°,  the  jet  is 
slender  anu  parabolic.  It  becomes  thicker  as  A  decreases,  and  reaches  a 
a  limiting  form  at  A  ■  0.  Then  there  is  a  stagnation  point  at  the  crest, 
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•"•iv-re  the  surface  makes  a  120®  a:igl-  with  itself.  This  angle  is 
dieted  by  the  same  argument  Stokes  used  in  his  study  of  water  waves. 

The  problem  is  formulated  as  as.  integro- different ial  equation  for 
the  two  free  surfaces  of  the  jet.  This  equation  is  discreticed  to  yield 
a  set  of  nonlinear  equations  which  are  solved  numerically  by  Newton's 
method.  In  addition  asymptotic  results  for  X  large  are  obtained 
analytically.  Graphs  of  the  results  are  presented. 


67.  "Fluid  dynamics  and  the  tie i r.rmnnn  equation"  by  Russel  E. 

Caflisch,  Stud.  Stat.  Mech.,  to  appear. 

The  fluid  dynamic  equations  and  the  Boltzmann  equation  provide  al¬ 
ternative  descriptions  for  the  evolution  of  a  gas.  In  the  limit  of 
small  mean  free  path,  the  Euler  or  Navier-Stokes  equations  can  be  de¬ 
rived  from  the  Boltzmann  equation  through  the  Hilbert  ot  Chapman- Enskog 
expansions.  These  expansions  are  singular  near  boundaries,  shocks,  and 
general  initial  data  and  special  boundary,  shock,  and  initial  layers  are 
needed  to  complete  the  solution.  In  the  outer  region  away  from  these 
layers,  the  expansions  are  shown  to  be  valid,  under  certain  conditions. 
They  approximate  a  solution  of  the  Boltzmann  equation  for  a  reasonable 
time  period,  and  they  provide  the  dominant  term  in  the  long  time 
asymptotics  of  the  solution  of  the  near  linear  problem.  Shock  wave  solu 
tions  of  the  Boltzmann  equation  are  also  discussed.  The  Boltzmann  solu¬ 
tion  for  a  weak  shock  is  constructed  and  is  shown  to  agree  with  the 
Navier-Stokes  profile  to  leading  order. 
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68.  "Theory  of  periodic  structures  in  lipid  bilayer  re-.arriaea" 
by  Meira  S.  Falkovitz,  Michael  Soul,  iiarry  L.  Frisch,  and 

Harden  M.  McConnell,  Proc,  Nat,  Acad.  Sci.,  submitted. 

An  approximate,  new  model  for  the  structure  of  the  periodic,  un¬ 
dulated  P  ,  phase  of  phosphatidylcholine  bilayers  is  proposed.  The  pro- 
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perties  of  this  phase  are  deduced  by  minimizing  a  Landau-de  Gennes 
expression  for  the  bilayer  free  energy  when  this  free  energy  contains  a 
term  favoring  a  spontaneous  curvature  of  the  membrane.  The  theoretical 
calculation  leads  to  a  model  for  the  Pg ,  phase  of  phosphatidylcholine 
bilayers  having  a  number  of  novel  physical  properties,  including 
periodic  variations  in  membrane  "fluidity." 

69.  "The  fluid-dynamic  limit  of  a  model  Boltzmann  equation  in  the 
presence  of  a  shock"  by  Russel  E.  Caflisch,  Institut  National  de 
Recherce  en  Informatique  et  en  Automatique,  INRIA  Report  No.  81,  June 
1981,  1-34. 

The  Broadwell  equation  is  the  simplest  model  of  the  Boltzmann  equa¬ 
tion  of  kinetic  theory  for  which  the  corresponding  model  fluid  dynamic 
equations  are  non- trivial.  For  this  equation  there  is  a  complete 
existence  theory  for  the  initial  value  problem.  Here  we  show  formally 
that  if  the  model  fluid  dynamic  equations  can  be  solved,  then  the  Broad- 
well  solution  asymptotically  converges  to  the  fluid  dynamic  solution  as 
the  mean  free  path  goes  to  zero.  This  limit  is  valid  even  if  there  is  a 
shock  in  the  fluid  flow,  although  there  is  a  thin  shock  layer  in  which 


the  convergence  Joes  not  held.  Arbitrary  smooth  ir.it  i:ii  u  »ta  .  ■  ,:i  lowed, 
which  leads  to  a  sliort  initial  layer  of  non-convergence,  hot  the  initial 
and  shock  layers  Jo  not  interact  due  to  the  assumed  initial  st.c cthness . 

61.  "Folds  in  solutions  of  two  parameter  systems  and  their  calcula¬ 
tion:  Part  1"  by  Allan  Jepson  and  Alastair  Spence,  Stanford  Univ.  Numer. 

Anal.  Report,  in  press. 

This  paper  treats  the  existence  of  paths  of  turning  points  in  solu¬ 
tions  of  nonlinear  systems  having  two  parameters.  It  is  well  known  that 
these  paths  are  solutions  of  a  particular  extended  system  of  nonlinear 
equations.  In  this  paper  both  regular  and  simple  turning  points  in  the 
extended  system  are  related  to  the  local  geometry  of  the  solution  surf- 
face  of  the  original  nonlinear  system.  A  description  is  given  of 
numerical  methods  both  for  solving  the  extended  system  and  for  calculating 
parameters  determining  the  local  geometry  of  the  solution  surface. 
Applications  to  perturbed  bifurcation,  to  the  formation  of  isolas,  and  to 
the  calculation  of  the  multiplicity  of  solutions  are  also  discussed. 


IV.  KGR:  PROGRESS. 


1.  "Optimal  catalyst  distribution  in  a  membrane"  (Joseph  B.  Keller, 
Meira  Falkovitz  and  Harry  Frisch),  sit  is  shown  how  to  distribute  a 
fixed  amount  of  catalyst  across  the  thickness  of  a  plane  membrane  to 
maximize  the  rate  of  reaction  of  a  substance  diffusin  through  the  mem¬ 
brane  and  reacting  within  it.  The  result  is  that  all  the  catalyst  should 
be  located  at  the  front  or  high  concentration  face  of  the  membrane.  It 
is  also  shown  how  to  distribute  the  catalyst  to  minimize  the  flux  of  un¬ 
reacted  substance  out  of  the  membrane.  In  this  case  the  catalyst  should 
be  uniformly  distributed  in  a  layer  which  is  symmetric  about  the  midplane 
of  the  membrane.  The  thickness  of  the  layer  is  determined.  The  analysis 
employs  the  methods  of  the  calculus  of  variations. 

2.  "Crawling  of  worms"  (Joseph  B.  Keller  and  Meira  Falkovitz). 

§The  mechanics  of  a  worm  crawling  along  a  flat  surface  is  analyzed.  The 
external  forces  of  friction  and  gravity,  and  the  internal  tension,  are 
taken  into  account.  An  equation  of  motion  is  formulated,  and  solutions 
are  sought  in  which  both  the  tension  and  the  linear  density  are  required 
to  lie  between  prescribed  bounds.  Periodic  travelling  wave  solutions  on 
worms  of  finite  length  and  pulse-like  travelling  waves  on  worms  of  in¬ 
finite  length  arc  treated.  The  maximum  crawling  velocity  is  determined, 
as  well  as  the  wave  form  which  achieves  it.  The  results  are  compared 
with  experimental  observations. 
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5.  "A  relation  between  the  ha  es  and  maximum  li.*  el  ihood  methods" 
(Joseph  B.  Keller).  §The  Bayes  transformation  B^_,  based  upon  an  observed 
value  y,  converts  a  prior  distribution  Pq(^)  of  a  parameter  0  into  a 
posterior  distribution  B  Pq(9).  J.  J.  Higgins  [Bayesian  inference  and 
the  optimality  of  maximum  likelihood  estimation.  Inter.  Stat.  Rev.,  £5 
(1977),  9-11]  has  shown  that  B^Pg(e)  converges,  as  n  increases,  to 
[9  -  9j(y)]  where  9^(y)  is  the  unique  maximum  likelihood  estimate  of  0 
based  upon  y.  We  show  that  when  there  are  several  maximum  likelihood 
points  in  the  support  of  p^S),  B^p^(9)  converges  to  a  sum  of  delta  func¬ 
tions  located  at  these  points.  Each  is  multiplied  by  a  certain  coeffi¬ 
cient  which  we  determine.  This  suggests  a  way  of  extending  the  maximum 
likelihood  method  to  the  case  of  several  maximum  likelihood  points. 


4.  "Weakly  nonlinear  geometrical  optics"  (John  K.  Hunter  and 
Joseph  B.  Keller),  §A  theory  of  weakly  nonlinear  solutions  of  systems  of 
hyperbolic  partial  differential  equations  is  contructed.  The  theory 
leads  to  propagation  of  waves  along  the  rays  of  the  corresponding 
linearized  theory.  However,  the  amplitude  satisfies  a  nonlinear  equation 
along  these  rays.  This  leads  to  waveform  distortion  and  shock  formation. 
The  theory  can  be  used  to  describe  weak  shock  waves  in  any  number  of 
dimensions.  Its  relation  to  the  theory  of  G.  B.  Whitham  [Linear  and  Non¬ 


linear  Waves,  Wiley,  1974]  and  of  L,  D.  Landau  is  demonstrated. 


"Helmhol t z 
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r-'sjn.itor  modes  and  scattering"  (Marian  3.  Sexton  and 
Joseph  B.  kellerj.  .  .-.n  integral  equation  is  formu  la  ted  for  the  problem 
of  scattering  of  r.n  acoustic  wave  by  a  Helmholtz  resonator.  The  equation 
is  solved  asymptotically  when  the  hzle  is  small  compared  to  all  other 
lengths  in  the  problem.  Scattering  resonances  are  found  to  occur  at  cer¬ 
tain  complex  frequencies  of  t he  incident  wave.  These  are  shown  to  be  the 
frequencies  of  the  modes  of  the  resonator.  These  modes  are  determined 
asymptotically.  The  corresponding  results  are  also  developed  for  the 
two-dimensional  problem  which  represents  water  waves  incident  upon  the 
narrow  mouth  of  a  harbor. 


6.  "Water  wave  production  by  an  oscillating  body"  (Joseph  B.  Keller 
and  Si-Xiong  Chen).  §Asymptotic  methods  are  employed  to  calculate  the 
waves  generated  by  an  oscillating  body  floating  on  the  surface  of  a 
fluid.  The  method  of  geometrical  optics  is  used  to  describe  the  far 
field,  and  boundary  layer  methods  are  used  to  describe  the  near  field. 

The  results  generalize  those  of  other  authors. 

7.  "Reflection,  transmission,  absorption  and  scattering  of  acoustic 
waves  by  rough  surfaces"  (John  G.  Watson  and  Joseph  B.  Keller).  S-The 
first  and  second  moments,  i.e.  the  coherent  field  and  the  two-point  two- 
time  correlation  function,  are  calculated  for  the  acoustic  fields 
scattered  from  various  rough  surfaces.  For  each  surface  they  yield  the 
reflection,  transmission,  absorption  and  differential  scattering  coeffi¬ 
cients,  as  well  as  an  equivalent  boundary  condition  for  the  coherent 
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;  i  l J  .  normalized  coefficients  arc  constructed  to  eliminate  divergences 
at  grazing  incidence.  The  results  are  specialized  to  surfaces  which  are 
statistically  homogeneous  in  both  space  and  time,  to  surfaces  which  are 
not  moving,  to  surfaces  which  are  simply  or  multiply  periodic,  and  to 
surfaces  consisting  of  randomly  placed  bosses  on  a  smooth  surface.  The 
surfaces  considered  are  slightly  rough,  moving,  soft  or  hard  boundaries, 
and  flat  surfaces  with  random  admittance  or  impedance.  The  analysis  is 
based  on  the  regular  perturbation  method.  Comparisons  with  previous  re¬ 
sults  are  made. 

8.  "The  prouuci-expansion  technique"  (John  G.  Watson).  §We  present 
a  method,  the  product-expansion  technique,  which  is  an  alternative  to  the 
smoothing  method.  Given  equations  for  a  random  process,  the  technique 
produces  an  approximate  solution  for  the  expected  value  in  the  form  of  the 
product  of  averages  of  operator  products.  An  application  to  rough  surface 
acoustic  scattering  is  given. 

9.  "Branches  of  periodic  solutions  and  their  computation" 

(Allan  D.  Jepson,  Esabius  Doedel  and  Herbert  B.  Keller).  §Branches  of 
periodic  solutions  to  autonomous  systems  of  differential  equations  with  a 
free  parameter  are  studied.  A  bifurcation  analysis  is  presented  with  a 
view  towards  numerical  methods.  A  computational  package  for  following 
these  branches,  detecting  bifurcation  points,  switching  branches  at 
simple  bifurcation  points,  and  switching  branches  at  Hopf  bifurcation 
points  is  tested  on  several  examples. 


10.  "Asymptotic  '’oundary  cend  L  t  ions  tor  ordinary  differential 
equations.  Parts  I  and  II"  (Allan  D.  Jepson  and  Herbert  5.  Keller). 

3 The  numerical  solution  of  no  point  boundary  value  problems  on  semi¬ 
infinite  intervals  is  often  obi lino !  by  truncating  the  interval  at  some 

finite  point.  In  this  paper,  we  determine  a  hierarchy  of  increasingly 
accurate  boundary  conditions  for  the  truncated  interval  problem.  Both 
linear  and  nonlinear  problems  are  considered.  Numerical  techniques  for 
error  estimation  and  the  determination  of  an  appropriate  truncation  point 

are  discussed.  Several  numerical  examples  are  .ivch. 

11.  "Fredholm  theories  for  two  point  boundary  value  problems  on 
semi-infinite  intervals"  (Allan  D.  Jepson).  §A  simple  method  to  deter¬ 
mine  the  Fredholm  properties  of  two  point  boundary  value  problems  on 
semi-infinite  intervals  is  presented.  The  general  problem  is  reduced  to 
an  equivalent  regular  problem  on  a  finite  interval.  The  Fredholm  pro¬ 
perties  then  follow  from  the  well  known  theory  for  regular  problems  and 
from  asymptotic  estimates  on  particular  solutions  of  the  singular 
problem. 

12.  "Paths  of  simple  bifurcation  points  and  their  computation" 
(Allan  D.  Jepson  and  Alastair  Spence).  §N'umerical  methods  for  following 
simple  bifurcation  points  in  nonlinear  systems  with  three  parameters  are 
presented.  Critical  points  in  these  branches  are  studied. 


14.  "Bifurcation  and  stabiiit.  in  constrained  variational  pro-..  ■' 
(John  H.  Maddocks).  §The  main  result  presented  is  a  theorem  determining 
when  a  quadratic  form  subject  to  linear  constraints  is  positive-definite 
Such  quadratic  forms  arise  as  tests  for  stability  in  general  constrained 
variational  problems.  Connections  with  bifurcation  and  exchange  of 
stability  are  also  described. 

15.  "Stability  of  the  elastica"  (John  H.  Maddocks).  §The  elastica 
is  a  standard  model  for  the  buckling  of  an  elastic  rod.  This  paper  de¬ 
termines  stability  properties  of  various  planar  equilibria  of  a  non- 
linearly  elastic  rod,  both  to  variations  in  the  plane  and  out  of  the 
plane.  Previously  undiscovered  second  bifurcations  with  associated  loss 
of  stability  are  described. 


